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Abstract—The problem of high-frequency diffraction by 
elongated bodies is discussed in this work. The classical 
approach is shown to require too high frequencies for its 
applicability. Attempts to improve the approximating properties 
of the asymptotic methods are discussed, focusing on the case of 
strongly elongated bodies examined in detail. The asymptotics 
are governed by the elongation parameter, the ratio of the 
longitudinal wave dimensions of the body to its cross-section. 
The cases of axial incidence as that of incidence at a grazing 
angle to the axis are considered. Both the asymptotics of the 
induced currents on the surface and of the far field amplitude 
are developed. Comparison with numerical results for a set 
of test problems shows that the leading terms of the new 
asymptotics provide good approximation in a uniform manner 
with respect to the rate of elongation. 


Index terms—High-frequency asymptotics, diffraction, elon- 
gated bodies. 


I. INTRODUCTION 


Methods of high-frequency diffraction remain an important 
tool for the analysis of wave phenomena in many applications. 
Classical results of Fock, Keller and others are restricted to 
geometries in which there is a single large parameter kp, which 
measures the characteristic size in wave lengths. The condition 
of validity of these asymptotic expansions requires that all 
other quantities describing the problem should not compete 
with this large parameter. In some cases, this imposes condi- 
tions that are too restrictive on the frequency, by requiring it to 
be very high. One of such cases is the diffraction by elongated 
bodies. 

Interest in high-frequency diffraction developed many years 
ago to fulfill the needs of numerous applications in radio- 
engineering. Since computers were not powerful at that time, 
alternatives to purely numerical analysis were highly appre- 
ciated. One of the important advances made on the basis of 
canonical problems was formulated in terms of the localization 
principle [1], which made it possible to find solution of 
diffraction problem for a specific convex body and deduce 
approximate formulas for the field distribution on its surface. 
These formulas are also applicable to any other body with 
the same values of the principal curvatures at a given point. 
Thus it is possible to decompose the domain occupied by 


the scatterer into smaller parts and to approximate each part 
with a canonical problem, whose solution leads to the local 
asymptotic representations valid in each of the subdomains. 
The matching of the local expansions is a nontrivial problem, 
albeit a solvable one as a rule. 

Application of the localization principle enabled Fock to 
find general formulas for the field in the penumbra region by 
working with the paraboloid of revolution as the particular 
body [2]. However, the same formulas can be obtained more 
easily with the help of the parabolic equation method. We will 
remind the outline of these derivations in the following section, 
since many of the ideas that form the basis of these derivations 
remain unchanged in the modern approach to solving the 
problems of diffraction by elongated bodies. 


II. FOCK ASYMPTOTICS 
A. The main assumptions 


Let us consider the wave field in the near vicinity of the 
light-shadow boundary on the surface of a convex body. This 
so-called Fock domain appears to be the cradle of creeping 
waves that propagate to the shadow region of the boundary, 
and of the Fresnel transition field in the penumbra region (see 
Fig. 1). 

Our first step is to introduce a coordinate system in which 
the surface of the body coincides with one of the coordi- 
nates surfaces. Usually such a coordinate system is formed 
by adding the normal to the semi-geodesic coordinates on 
the surface. That is, we consider a set of geodesics whose 
directions are defined by the rays of the incident field, and 
measure the first coordinate s along these geodesics starting 
from a reference line, which is usually the light-shadow 
boundary on the surface. The second surface coordinate a just 
parameterizes the geodesics from this set. In this approach, 
how the normal coordinate is defined is not important and 
it can be taken as the geometrical distance to the surface 
for the case of a general body. However, when dealing with 
elongated bodies we will return to this question of the choice 
of the normal coordinate. Once the coordinate system is 
chosen all the equations and the boundary conditions of the 
problem should be rewritten in this coordinate system. The 
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Fig. 1. Diffraction by a convex body. Fock domain — red, Fresnel penumbra 
— blue, creeping waves — green. 


characteristics of the surface and the incident field, such as 
radius of curvature p of the geodesics, their divergence h, 
radius of transverse curvature p, and the torsion appear in the 
above equations. 

The next step in this procedure is to choose the scales of the 
coordinates. One should scale the coordinates in a way such 
that the derivatives of the field by the scaled coordinates can 
be considered as quantities that are on the order of unity. In 
our case, the field depends on s in the form of a dominant 
factor e’*$, multiplied by a slower varying function called 
the attenuation function. To define the scales one can proceed 
in one of two ways. One possibility is to search for these 
scales in the solution of a canonical problem, which can be 
solved exactly. The other possibility is to look at the terms 
in the equations and to attempt to choose terms that should 
be considered as the principal ones. We will illustrate this 
second approach by considering the diffraction by a cylinder 
as an example. For this 2D problem, Maxwell’s equations are 
reduced to the Helmholtz equation for scalar waves. 

Let S be the convex boundary of a body illuminated by a 
wave field given by its ray expansion. Figure 1 presents the 
two-dimensional cross-section of the field of rays. On the line 
C on the surface S' (see Fig. 1) the rays of the incident field are 
tangent to S. Such rays are called limiting rays. The surface is 
illuminated on one side of C, that is it is reached by incident 
rays and each incident ray is reflected from S. No geometrical 
rays reach the surface on the other side. 

As mentioned above, our next step is to we introduce the 
representation 

U = exp(iks)u(s, n), (1) 


and assume that the new unknown function u varies more 
slowly with the s coordinate than does the exponential multi- 
plier. We substitute the above representation into the equation 
and rewrite it in terms of the coordinates (s,n) as 








Ou Ou n pf. Ou Oru n Ou 
21k iku 4 
ðs ðs? p+nop ðn On? p On? 
n? 3u lôu nôu n n? 
! + 2k?-ut+tk?—u=0. (2 
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Next, we introduce the stretched coordinates 
ks, kfn 


where the exponents a and ĝ are unknowns as yet. We 
assume only that a < 1 (otherwise u varies with s faster than 
the exponential factor in (1)). This assumption enables us to 
exclude the term 0?u/0s? from principal terms and to consider 
it as a correction to the term 27kOu/0Os. Further, to be able to 
fix the boundary condition at n = 0 and radiation condition 
at n = +00 the second-order derivative w.r.t. n should be 
included in the principal-order terms. This implies that 6 > 0 
and the terms with Ou/On are smaller than 0?u/On?. Finally, 
the n? terms are smaller than similar terms with n in the first 
power. In view of all this the following three terms may be 
considered to be the principal terms in equation (2) : 
2 

Be aye 

On? Os 
One should verify that all the terms in (3) are principal- 
order terms, because otherwise the boundary value problems 
involving n would have no solution. Preserving all three terms 
yields the well known Leontovich parabolic equation : 


ðu Ou n 
33+ pnz | 2k a (4) 
Simultaneously, we have determined the size of the Fock 
domain. Parabolic equation (4) fixes the scales a = 1/3 and 
B = 2/3. Thus, Fock domain is as small as k~!/3 along the 
surface, and as small as k~2/3 in the direction of the normal. 
For the more complex case of 3D problem, the parabolic 
equation remains the same, since within such a small domain 
any general surface is almost cylindrical. 


and 2k u, (3) 


2ik 





B. Solution of the parabolic equation 

In the Fock domain, we can usually replace p(s, a) by its 
value po(a@) on the light-shadow boundary, since the Fock 
domain is small and p(s, œ) ~ po(a). Then the coefficients of 
the parabolic equation become constant, and it can be solved 
by means of Fourier transform to get : 


+00 
“= f aly, Ode, (5) 


where we have introduced the stretched coordinates 


s n k oi 
o = m—, v = 2m? —, m = (==) . (6) 
po Po 2 
For u, we obtain the Airy equation 
où ` 
za + (¥-6)u=0. (7) 


We use the Airy functions v and w, in Fock notation. 
The function wı is chosen because it satisfies the radiation 
condition representing only the waves which propagate away 
from the surface. The solution can be written as an integral 


+oo 


u= / ei e(C) (o - v) + R(d)wi(e — v) )ac. (8) 


— 00 
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The part of the solution which contains the Airy function 
v(¢ — v) corresponds to the incident field and the part of the 
solution which contains the Airy function w;(¢—v) represents 
the secondary field. The coefficient R is determined from the 
boundary condition. Thus from the Neumann condition, for 
the TM polarization we find 

o(¢) 
w (C)? 
where dot denotes the derivative. For the TE polarization, by 
imposing the Dirichlet condition we find that 





R™(¢) =- (9) 





TEA vd) 
i (6) = wi(¢) 


The amplitude c(¢) remians undetermined as yet, and it 
enables us to match the local asymptotics in the Fock domain 
to the incident field. For that we consider the incident wave and 
represent it in the stretched coordinates (ø, v) of the boundary 
layer. For example, for the case of an incident plane wave we 


have i 
ut = u'(C) exp {i (ov — 57°) } . 


By inserting this expression in (8), we obtain the integral 
equation for the amplitude c(¢), which reads 


(10) 


+00 


f P A E O fi (ov £ 50) \ . 


a a1 
Applying the Fourier transform we find 


; 1 
c(¢) = Ei 
If o is large, one can compute the integrals by using 
the residue theorem. Then, the contribution of every residue 
matches to the creeping waves which propagate along the 
surface of the body. 


C. Induced currents: leading- and higher-order terms 


If we take the observation point in formula (8) on the 
surface, i.e., if we let v = 0, we obtain the well known 
expression for the field associated with the TE-polarized wave 


giei 


Similarly, for the induced current of the TM wave, we get 


J= aE to AG 

OF th 
Although (12) and (13) only contain the leading order 
terms of the asymptotics, and the asymptotic decomposition is 
carried out by using the inverse powers of m, i.e., fractional 
powers of kp, they yield a rather accurate approximation 
for the field even when the frequency is not too high. We 


can verify this from Fig. 2 of [3] (the original is in [4], 
where the currents induced on a circular cylinder, whose 


E = U’(C 





(12) 





(13) 


diameter is equal to the wavelength, have been presented. It 
is worthwhile to point out that in the diffraction problem the 
current is formed by the interference of two waves. One travels 
in the clockwise direction from the light-shadow boundary, 
while the other originates from the light-shadow boundary on 
the opposite side of the cylinder and travels in the counter- 
clockwise direction. Both of these waves are attenuated, but 
as the diameter of the cylinder is not large, the wave which 
encircles the shadowed side, does not travel a large distance, 
and its contribution remains noticeable in the form of the 
oscillations of the field amplitude. 


For an arbitrary (non-cylindrical) 3D surface, the derivations 
of the principal-order term are exactly as presented above. 
The second surface coordinate œ is not stretched in the Fock 
domain and derivatives with respect to this coordinate are 
not involved in the higher-order operator. Because of this, the 
leading-order term of the asymptotics is not affected by the 
transverse curvature of the body. However the induced currents 
are different if the body is a sphere, and the agreement between 
the real currents and the asymptotics is not so good as it is 
for the case of the cylinder. A comparison with the numerical 
results presented in [5] defines the domain of applicability 
where the results from the Fock formulas suffer from less 
than 5% error, provided kp > 20. Attempts to improve the 
situation entail the computation of the next-order terms of 
the asymptotics, and we mention here only three papers that 
have attempted to do this. The first one is by Hong [6], 
who has derived the second-order term of the asymptotics, 
both for the fields in the Fock domain and for creeping 
waves in the deep shadow. He has considered an arbitrary 
(with some restrictions) smooth convex surface with ideal 
boundary conditions. For the impedance boundary conditions, 
the second-order term in the asymptotics of creeping waves 
was found in [7]. It was possible to find not only the second- 
order term for the case of a sphere, but the entire series 
representation [8]. However, the latter result differs by a 
factor of 2 from the one by Hong. Inspite of the possibility 
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Fig. 2. Currents on a perfectly conducting circular cylinder with the diameter 
equal to the wavelength. Solid line — exact, dashed — Fock asymptotics (13). 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


to use the second- and higher-order terms, the lower bound 
of frequency domain where the asymptotic formulas yield 
sufficiently accurate representations for the currents on the 
sphere, have remained almost unchanged. This leads us to 
conclude that the higher-order terms do not enlarge the domain 
of applicability of the first-order representation. 


III. MODERATELY AND STRONGLY ELONGATED BODIES 


For elongated objects, such as spheroids, the lower fre- 
quency bound of the domain where Fock formulas can be 
used shifts to kb > 100 and larger. We can deduce this 
by comparing the principal-order term of the asymptotics 
with numerical results given in [9]. Figure 3 presents the 
induced currents on the spheroid with semiaxes a = 0.5m 
and b = 1.25m, at frequencies 1, 2 and 4 GHz and the 
asymptotics. We use the stretched coordinate o in which 
the “size” of the spheroid changes with frequency, but the 
asymptotics remain unchanged. Figure 3 shows that the Fock 
formula underestimates the values for the currents. This is 
consistent with the observation that the edges or surfaces 
with high transverse curvature decrease the attenuation of 
the creeping waves. That is, large transverse curvature or 
sharp edges promote the propagation of waves, as exemplified 
by the following cases. While studying diffraction by disks 
Senior [10] observed, in 1969, that a wave propagates with 
low attenuation along the rim of the disk. The Sommerfeld 
wave [11] propagates with only logarithmic attenuation along 
conducting wires. Finally, a source on a cylinder excites a 
wave that propagates along the generatrix of the cylinder, and 
only decreases as the inverse of the square root of the distance 
[12]. 

The impact of transverse curvature on the propagation 
constant of creeping waves, and more specifically, the case 
of creeping wave on elongated bodies, has been analyzed in 
[13] (see also [14]). In the spirit of the work of Engineer 
et al. [15], who have studied the problem of diffraction by 
2D slender bodies, it was assumed that besides the usual 
asymptotic parameter kp, the geometrical characteristics of 
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Fig. 3. Induced current amplitude on the spheroid with semiaxes b = 1.25 m, 
a = 0.5 m at frequencies 1, 2 and 4 GHz (curves 1,2,3) and the asymptotic 
approximation by Fock (dotted curve). 


the body form another large parameter A = p/p;, with 
which the asymptotic parameter may compete. Analysis of 
the scales shows that reorganization of the usual creeping 
waves asymptotics as well as the asymptotics in Fock domain 
takes place when A reaches the order of m, i.e., becomes 
as large as (kp)'/3. This case was referred to in [13] as the 
case of moderately elongated body. For this case, the leading 
term of the asymptotics of creeping waves starts to depend 
on the transverse curvature, and transverse curvature on the 
order of (kp)'/3/p changes the creeping waves asymptotics, 
although these changes are not too significant. The structure 
of the asymptotic expansion remains relatively unchanged, 
though some terms in the recurrent system of the boundary- 
value problems move from the second-order corrections to the 
first-order. This results in an aditional multiplier defining the 
amplitude of creeping waves. This result can be simulated by 
introducing an effective impedance : 
a 

E 2kp: 
which can evidently be used in the Fock domain. The asymp- 
totics of induced currents on a moderately elongated body 
improves the approximation, as may be seen from Fig. 4. 
However, the currents still remain underestimated. 

The strongly elongated body was defined in [13] as the body 
in which the transverse curvature is on the order of k?/3 p713, 
implying that A ~ m?. The asymptotics of creeping waves on 
such a body changes significantly and leads to the equation 
below, which is more complex (see [13]) than the Airy 
equation (7) 





(14) 


faba 3 OU 
Ss —¢€)a=02 15 
Ov? n tte (ja : (15) 
where 
KZ on? Pe (16) 
p 


For the case of a strongly elongated body, the parameter « is on 
the order of unity and is the main parameter which describes 
the effect of the transverse curvature. The formulas reduce 
to the effective impedance approach of moderately elongated 
case, if we let the parameter « to be small. Equation (15) is 
a biconfluent Heun equation [16] which has 4 singular points, 
three of which coalesce at infinity, while the last one goes to 
zero. The singularity of v = —« in (15) corresponds to the 
points on the axis of transverse curvature and it is convenient 
to introduce a “shifted” normal v’ = v + «k, especially for 
the case of a body of revolution. Then the surface equation 
becomes v’ = «. In order to be able to write the asymptotics of 
creeping waves on a strongly elongated body, it is necessary to 
impose some restrictions on the variation of both the curvature 
and the transverse curvature; specifically it is necessary to 
assume that both p and p; vary very slowly. 

Attempts to generalize the approach to Fock domain (almost 
the only published result we can locate is [17]) were neither 
mathematically rigorous nor satisfactory from the point of 
view of providing a good approximation for the currents. 
Although the effect was apparently described correctly at 
the qualitative level, agreement with the results of numerical 
computation was less than satisfactory. 
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Fig. 4. Moderately elongated body approximation for diffraction by the spheroid with a = 0.5m, b = 1.25m at 1 GHz (left) and 2 GHz (right). Numerical 
results — solid, moderately elongated approximation — dashed, Fock approximation — dotted. 


It was understood later that the restrictions imposed on 
the curvature variation were inconsistent with the geometry. 
Indeed for a very elongated body, such as a prolate spheroid 
with A = O(m?), which is illuminated along its axis by 
an incident wave, the entire object is not located in the 
deep shadow region, but in the light-shadow transition zone. 
Therefore it is not correct to replace it with an object with 
constant p and pz. 

All the above-mentioned facts led to the conclusion, which 
was recently rediscovered in [18], that a different technique 
is neded to handle the case of strongly elongated body. To 
address this problem, a new technique was introduced in 2009 
[19], [20] which is currently being developed. The rest of the 
paper describes this new approach. 


IV. DIFFRACTION BY A STRONGLY ELONGATED BODY 


A. Boundary layer, scales and coordinates 


Let the elongated body have the symmetry of revolution 
and let its main radii of curvature at the widest cross-section 
be p and p+. We assume that the body is spheroid-like, i.e., 
its surface deviates from the surface of the spheroid by an 
asymptotically small distance. The radii of curvature define 
the semiaxes of this spheroid as 


b= ppt. 


The assumption that the body be strongly elongated in terms 
of the semiaxes can be expressed as 


a= pt, (17) 


a ae O(1). (18) 


We refer to the parameter x as the elongation parameter. It is 
related to the parameter «, introduced earlier, via the formula 


x= 9-1/2 3/2, 





Fig. 5. Geometry of the problem. 


We use the spheroidal coordinates (€,7,y) related to the 
cylindrical (r, z,) ones, where the z axis is the axis of the 
body (see Fig. 5), via the formulas 


2 py. TS py -1y1-n?, 


where p = yb? — a? is the semi-focal distance. 

In order to choose the scales, we initially consider the case 
of a plane wave incident along the axis of the body, and make 
some generalizations later. If the body were not elongated, 
we would use Fock asymptotics described in Section II. The 
size of Fock domain is on the order (kp)~'/° in the direction 
of the wave incidence, while it is on the order (kp)~?/% 
along the normal direction. We see that if we maintin these 
orders, then the entire strongly-elongated body would be inside 
the Fock domain, implying that we should not stretch the 
angular spheroidal coordinate 7, but should do so with the 
radial coordinate €. It is convenient to replace € with another 
coordinate, which we denote by 7, such that its value measured 
from the axis would be on the order of unity on the surface. 
We introduce such a coordinate by using the formula 


(19) 


T = 2x7 'kb(é — 1). (20) 


So that 7 = 0 on the axis while it equals 1 on the surface. 
The coordinates (n, T) serve as the coordinates of the boundary 
layer, and they replace the usual o and v coordinates of the 
asymptotics of Fock. 
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Keeping in mind that a < b, we can simplify the formulas 
for the boundary-layer coordinates to 


r=aV/1—17y7, 
2 


z =bņ+ z? —1)n. 
This change does not influence the asymptotic procedure 
which pertains to the leading order, although may affect the 
higher-order corrections, which we ignore however. 

For an ordinary (not elongated) body, the use of the 
parabolic equation method suggests that the factor e’** be 
separated. Computing the arc-length s on the spheroid is not 
so easy; however, we can use a simpler factor e’”?” or etk’ if 
we again exploit the fact that a < b. These factors agree with 
e’*s up to the terms on the order of unity, and we can choose 
any of them. The formulas become a little more compact if 
we choose e’*?", which we simplify to 


exp (iton — 5x) ‘ 


(21) 


(22) 


Once we have chosen the coordinates and the factor men- 
tioned above, we can go ahead and generalize the problem, 
by allowing the wave to be incident at an arbitrary angle with 
respect to the axis. However, the angle should be small, such 
that the phase shift corresponding to the trajectory of the wave 
along the geodesics on the surface differ from the phase shift 
of the multiplier (22) by no more than a quantity on the order 
of unity, implying that the angle of incidence 7 should be such 
that 


B = Vkbd = O(1). (23) 


Dimensionless quantities x and £, that are on the order of 
unity, are the parameters of the asymptotic formulas, which 
we will derive below. 


B. Integral representation of the field 


We search for the solution of Maxwell’s equations in the 
form of a Fourier series in terms of the angle y. For the part 
of the electromagnetic field which depends on the angle in 
the form e’“”, we express all the components of electric and 
magnetic vectors via the Æ, and H, components. Below we 
denote the Fourier harmonics of these components as Ez and 
Hz. The functions Ey and Hy are solutions of the system 
of differential equations, which is cumbersome to handle, 
although it is standard in spheroidal coordinates. Separation 
of variables is possible only if £ = 0 [21], and it leads to a 
representation of the field that utilizes spheroidal functions. 

We now extract the dominant factor in (22), and neglect 
the lower-order terms. For the leading order in terms of the 
asymptotic parameter kb, we obtain the system of parabolic 
equations. It can be easily split into two independent equations 
for the new unknowns as follows : 


Fig +t Ey — ig 


P(T, n) = 2 Q(T, n) = 2 (24) 
These unknowns satisfy the equations 
Le- P; = 0, Le41Qe = 0, (25) 





= cat. 2 
la Goat get a N 


2 
ri ler = La n) ). (26) 





4 


The differential equations in (25) are considered together with 
the boundary conditions. For the case of a perfectly conducting 
surface, the boundary conditions require that the tangential 
components of the electric vector be zero on the surface, which 
leads to the conditions 


P(1,7) + Qe(1,) = 0, (27) 
AP(1,n) 1 OQc(1,7) 1 m 


We also set the radiation conditions at infinity for the sec- 
ondary field. 

Parabolic operator (26) is amenable to a separation of 
variables. Elementary solution of the equation L,U = 0 can 
be written in the form 


1 1—n\" 
Un = TV1 — 1? (; 4 1) Fanj2 (—ixT) , (29) 


where F, n/2(g) is a solution of the Whittaker equation [22] 


1 l-n? 
ra ( ae P) F=0. 








4 g Ag? 


We choose the Whittaker functions M,,,/2(—ixT) for the 
incident field in (29). These functions are regular at 7 = 
0, which results in solutions that are regular in the entire 
space. For the secondary field, we choose Whittaker functions 
Wa, n/2(—ixT) which correspond to solutions satisfying the 
radiation condition in terms of T. 

We combine P; and Qg of the elementary solutions (29) 
in the form of integrals over the separation parameter u. In 
these representations we include amplitude factors that are 
dependent on u. To define the path of integration, and the 
amplitude factors, we consider below the incident field first. 


C. Representation of the incident field 


An arbitrarily polarized plane electromagnetic wave can be 
represented as a sum of transverse electric (TE) and transverse 
magnetic (TM) waves. We define the transverse direction with 
respect to the plane of incidence and assume that it coincides 
with the X OZ plane in Cartesian coordinates (see Fig. 5). The 
plane of incidence is not defined for the axial case, but TE and 
TM waves can still be defined with respect to the plane XOZ. 

Let us consider the TE case first. We set a unit amplitude 
for the incident TE wave and express its E and H fields as 

E’ = exp (ikz cos 0 + ikx sin V) ey, (30) 

H’ = exp (ikz cos V + ika sin v) {— cos Vey, + sin vez}, 

(31) 
where ez, €y and e, are the unit vectors of Cartesian coordi- 
nates. Next, we represent the incident wave in the form of a 
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Fourier series as follows : 


E; = exp (ikz cos ù) fis (kr sin v) 


n=l 


+2 >D i”! J, (kr sin 8) cos(ny) \ (32) 


Jn (krsinv) . 
kr sin 0 amine): 
(33) 
Each of the harmonics of (32), (33) can be expanded in an 
asymptotic series in terms of inverse powers of kb. The leading 
order terms satisfy the parabolic equations (25). Therefore 
the leading order terms in the asymptotic representation of 
(32), (33) can be represented as a combination of elementary 
solutions (29) in which, as explained above, we should include 
the Whittaker function M, ,/2(—ixT). Expressions for the 
harmonics in (32) and (33) contain the derivatives of Bessel 
functions and the ratios of Besel functions to their agument. 
However both these expressions can be reduced to just the 
Bessel functions with the help of the formulas : 


2Jn(9) = In—1(9) — Int. (9); 


In(9) 
g 
from [22]. Therefore the problem of finding the leading-order 
representation of the incident wave as a combination of the 
elementary solutions (29) reduces to that of finding such a 

representation for the functions 


oo 
Hi =2exp (ikz cos 9) cos 0 5 iin 


2n = Jn—1(g) + Jn+1(9) 


exp (ikz cos 0) Jn (kr sin ¥). 


When substituting here the expressions for the coordinates r 
and z we replace the trigonometric functions of small J by 
their approximations : sind ~ V and cos) ~ 1 — $0?. By 
doing this, we obtain the equation 





1) Mn /2(—txT du = 
=V(n,7). G4 


Both the left-hand as well as the right-hand side of this 
equation satisfy the parabolic equation with operator Ln. 
Therefore, if we find the amplitude A„(u) for any fixed 7, 
it will satisfy (34) for all other values of 7. This important 
property simplifies the problem of finding the solution. First, 
we rewrite (34) in the form 


rca! (G4) 4 i 





p 
ma (Fey) MOM" 
My,n/2(-ixT dp = V(n,7), (35) 
where 
Myn/2(—i8?) 
Vx B 


We note that 7 is presented in (35) only as the product yr. 
Thus, since Q, does not depend on 7, it does not depend on 
x either. Further, since y7 and —p?, presented in (35), are 


An(H) = On (u). (36) 


symmetric, Q, does not depend on 8 either. This enables us 
to find Q, for any chosen x7 and 8. Equation (35) takes the 
simplest form if 8 and xT both > 0. Computing this limit 
and using the asymptotics 





Munp(g)~g"t?, 90, (37) 
1 sg\” 
Jalo) ~ (2) ) g>0 
we reduce equation (35) to 
=a 1 9) (n+1)/2 
—— ) 0, (p)dp = 1- ; 38 
E (u)du lA n’) (38) 


We remind that we have not yet chosen the path of inte- 
gration. So if we choose it now as the imaginary axis and 
set 

s—1 





s+1 


? 


equation (38) reduces to 


+ioo 
2 Vs n+l 

“#0, (udu = — 

ade = A (34) 


We identify the left-hand side of this equation with inverse 
Mellin transform and find 





—ioo 








1 n+1 n+1 
Q = T r — ; 
oi) = at (TE +n) (AE - a). 69 
Using (39) we can write the incident plane TE wave in the 
form : 
gi — _j&xplekbn — ixn/2) 





- TE 
? V1 = PB /XF gy 
x -osiota (—ixT)Mi1/2(i8°) 
+00 
+ 5 if cos(Ly) [Qe (it) Muy es 
l=1 


-Oera (it) My xs (—ixT) My ess (68?)] ) dt, (40) 


R bent 
l eo g 


M, e (i6”) 


e-1(—ixT)M it, 1 (i8? ) 


jcxp(ekbn = ixn/2) 


v1 -= n’ BSxT 


+00 
x O i sin(€y) [o (it)M,, « 
f=1 





H; = 


e-1 (7x7) 


wo (0B?) ] dt. aD 


The representations for the TM plane wave 


+0041 (it) Mj, e (—ixT)M, 
E’ = exp (ikz cos 0 + ikz sin V) {cos Ve, — sin vez}, 
H’ = exp (ikz cos 0 + ikz sin V) ey 


can be derived analogously. The electric field is given by the 
expression (41) and the magnetic field can be expressed as in 
(40) with the opposite sign. 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


D. Boundary conditions and the induced currents 


The representations for the secondary fields are similar to 
(40) and (41), the only change being the replacement of Whit- 
taker functions Mit n(—ixT) with the functions Wiz,(—ixT), 
multiplied by a function of the integration variable t. Note that 
the Whittaker functions W;t n(—ixT) are chosen because they 
satisfy the radiation condition. For the TE case we have, in 
accordance with (24) and (25), 


+00 


1 it 
[Ga 
1+ 


gs — _,oxp(ikon — ixn/2) 


o I= PBT 


x {BTEC (A 





+00 


+ 2 if cos( lip) [Qe (it) Rele) Wy e1 (—ix7) My, 2-1 (iB?) 
f=1 


v 


=e (it) Tet) Wp, 241 (—iXT) My, s (68?)] dt, 


(42) 
+00 


L—n\it 
| (aa) * 
1+ 





He = exp(ikbn — ixn/2) 


0 O JI- BY 


+00 
x X- if sin(Ey) [Qe (it) Re(t) Wy, 2-1 (ivr) My, 11 (iB?) 
f= 


Oesa (it) Te(t)W y eps (xT) M; t41 (1B”)| dt. (43) 
The functions Rọ(t) and T(t), introduced herein, play the role 
of the reflection coefficients in some sense, and can be found 
when the sum of the incident and secondary fields is substi- 
tuted in the boundary conditions (27), (28). Evidently, each 
of the harmonics satisfies the boundary conditions separately. 
For the TE case, we find (for details of the derivations see 
[23]) : 
i Mit,1/2(-tx) 











To = ae (44) 
° Wit 1/2(—ixX) 
Ye-1e41 1 Wa-1)/2 
Ty = l=1,2,... (45 
£ Ze Z Q ` 2, (45) 
Yei -1 1 
Re = a a = Z; Mern/aCes £=1,2,... (46) 
where = 
CG, = Etat Minern al) 
202(€ + 1)? Mie,ce-1)/2(87)’ 
Yn,m = Wi, g (—ix) Mit, m (—ix) 


+ Wit,» (—ix) Mit,m (—ix). 


Zn = Wae (—ix)W;y, 2 (—ix) 
+ Wix, e= (—tx)Wiy, a4 (—tx) 


t 


and W,, is the Wronskian of Whittaker functions [22] 


Wr = Mitn(g)Wit.n(9) Wa Mitn(g)Wit,n(g) 
— T+ 2n) 
T(1/2+n -— it) 





(47) 


For the TM case, the formulas are similar. For this case, 
E, is given by the right-hand side of (43) while H, is given 
by the right-hand side of (42), with the common multiplier 
(—1). The reflection coefficients in these two formulas are as 
follows: 

T) = _ Many( 8x) 
Wit 1/2(~¢xX) 


and the other coefficients Ty, 2 = 1,2, ... and coefficients Re 
are defined by formulas (45) and (46), respectively, but with 
Ce+1 replaced with —C¢e+1. 

The expressions given in (40) through (43) describe the 
electromagnetic fields in the boundary layer near the surface 
of the spheroid. In contrast to the classical Fock asymptotics, 
they do not have locality with respect to the transverse surface 
coordinate. Indeed, by restricting our analysis to the bodies 
of revolution we obtained the Fourier series in terms of the 
transverse angular coordinate ọ. In other respects, the structure 
of these set of formulas has some similarity with the classical 
asymptotics of the field in the Fock domain. For example, the 
Fourier integral in Ç in the classical asymptotics is replaced 
by another integral transform, and the Airy function v(¢ — v) 
which represents the incident plane wave in (8) is replaced 
with the Whittaker function M;t n(—ixT). The Airy function 
wi(¢ — v) which represents the secondary wave outgoing 
from the surface is replaced with the Whittaker function 
Wit n(—ixT). The reflection coefficients (9) and (10) became 
more compex and are given by the expressions in (44) through 
(46), and (48). 

The induced currents J = n x H|,—; can be found by using 
(40) through (43), with its principal component directed along 
7. After simple manipulations we can obtain 


(48) 








J =e" A(n, x, b, p), (49) 
where 


2 en ixn/2 


I= FV 


+4 , 
: fc — 2)" 5 if! sin(£y) S 
1l+7 = (€+1)!Z, 


£ 


A™®(n, x, b, p) = 








£ 
x fr aba is) My, ep (18°) W,, e (~x) 








£ 
H L(L+1)r (5 i) My, (i8)W en (ix) } dt 
in the TE case. 
It is also convenient to represent the special function AT? 
in terms of the Coulomb wave functions F and H~ [22], and 
to use the program developed in [24] for their computation. 
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Using the formulas 


My, ea (18°) 
7 2ietTt/4+7t/2P(£ + 2) 2 ( £) 











J¥(e/2+1+ir/2+1—it) O 2 
Witt (ix) 
a T(£/2 + it) x 
= pine /A+nt/2 Ht ( t ) 
A T(é/2—it) 3-1\ 72 
we obtain 
ee —8e-*xn/2 





“ROS it CO 
= i” sin(€y) x 
erie = 2, 














p 
F A + Fea (t, 5) HE 48) 
HELD ATA) CA 
(50) 
Similarly, for the TM case we can obtain 
8e~ixn/2 
A™ (n, x, b, p) = ——= X 
m/1— n? /XB 
a it F, (4 5) 
l-n 1 +0 
J G paap he 
2 Q2 
Fs (t, F) Wha (-t, 8) - Fea (t, D ALY), 
Hi, (-t,$) HT (-t, $) + Hike (-t, 3) H} (4 §)J 
(51) 


It is worth noting that one should consider the limit as 6 —> 
0 for the case of axial incidence. Only the Coulomb wave 
function f'-1 has a non-zero contribution in this case and the 


2 
special functions ATE and A™ simplify to 











ATE (n, x, 0, p) Z A(n, x) sin y, 
AL (n, x; 0, p) = A(n, x) cos y, 
where 
AG 4e*xn/2 G = n) t eTt/2 7 
úi VIXV1— 1 » Loe a cosh(zt) 
: Hy.(—t, %) dt 
H*,(-t,¥) Hy (-t, ¥)+H*, (-4, 3) HY (-+, 3) 


(52) 


The expressions in (50), (51) and (52) are more cumbersome 
in comparison to the classical Fock functions (13). However, 
with the help of the program from [24] their computation are 
relatively straightforward. The integrals converge quite rapidly 
and only a small interval contributes to the integral. As shown 
in [9] the function A reduces to the Fock function (13), when 
X 7 +00. 


E. Backward wave 


The total field in the boundary layer near the surface is 
the sum of the primary wave running in the positive direction 
of z together with waves that are formed when the primary 
wave encircles the ends of the spheroid. The asymptotics in 
(50) and (51) only describe the primary wave and are valid in 
the middle part of the spheroid. The reflected backward wave 
representation in the boundary layer can be written by noting 
that the backward wave runs in the opposite direction and that 
there is no incident backward wave. For the case of the axial 
incidence, considered in [25], the representation of the current 
corresponding to the backward wave is given by 


dg Se. BG a) sing, JE™ Se By) cosg, 











where 
T +ioo 2 
4eixn/2 È 1+7 it e77t/2 
TE (es 
JIxVv1i- J \t=n cosh(rt) 
Hy (~t, X) r(t) dt 
x - 3 
Ht, (—t, %) Ht (—t, $)+H7, (—t, $) HT (—t, %ž) 


(53) 


The special function B is defined by the expression similar 
to the one given in (52) in which we change the sign of 7 
and introduce an additional multiplier r(t). The latter may 
be interpreted as the reflection coefficient from the tip of the 
spheroid. To find this multiplier one needs to construct the 
field in the vicinity of the spheroid tip, that is for 7 ~ 1, and 
match the sum J + J, with that field. The difficulty is that 
there is no asymptotic parameter in the region near the tip of 
the spheroid. Indeed, the radius of curvature of the surface at 
the tip is pe = a?/b. Therefore the parameter kpe, which is 
equal to the elongation parameter x, is on the order of unity. 
However, there is a way to circumvent this difficulty, as shown 
in [25], which is based on the observation that a very elongated 
ellipse, such as a trajectory of a point mass in the central field 
of the gravity, is well approximated by a parabola. Hence, if 
we replace the surface of the spheroid by the surface of the 
paraboloid defined by the equation 


2 
r? =25-(b-2), 


we can use the exact solution of Fock [2] to approximate the 
field. Matching this Fock solution to the sum of the forward 
and backward waves allows the reflection coefficient to be 


found 
= je2tkh—-tx T(1/2 F it) na =2it (54) 
SPEER bee 





r(t) 


We note that the asymptotic parameter kb is presented in 
(54) and, rigorously speaking, the formula requires further 
asymptotic simplifications in the process of computing the 
integral in (53), as the contribution of the residues in the zeros 
of the denominator 


HY, (-4.5) #7 (-45) 
73 2 2 2 


+ Ht, ( 1) HY ( 1,5) =0. (55) 
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TABLE I 
TEST PROBLEMS PARAMETERS. 


























no. | f(GHz) | a (m) b (m) x kb 
1 1 0.5 1.25 4.19169 | 26.19806 
1 2 0.5 1.25 8.38338 | 52.39613 
2 1 0.5 1.76776695 | 2.96397 | 37.04966 
2 2 0.5 1.76776695 | 5.92795 | 74.09931 
3 l 0.3125 | 1.39754249 | 1.46452 | 29.29032 
3 2 0.3125 | 1.39754249 | 2.92903 | 58.58065 








Solutions of equation (55) lie in the lower complex half-plane 
of t (see [20] where solutions of (55) are analyzed) and the 
main contribution to (53) is due to the residue in the pole 
with the maximal imaginary part of t. However, finding the 
solutions of the dispersion equation (55) is not easy, and we 
compute the integral directly instead. Moreover, we expect that 
the formulas can also be applied in the case of values of kb, 
that are not too large, say kb ~ 3, when the approximation 
with one residue fails. 

To illustrate the approximating nature of the asymptotic 
formulas we present the Figs. 6 through 8 taken from the 
paper in [25]. These figures present the currents induced by 
the axially incident plane wave on the surface of perfectly 
conducting spheroids. The trigonometric multiplier of y is not 
considered. The parameters of the test examples are presented 
in Table I. The results computed by the Finite Element Method 
(FEM) are plotted in red and the asymptotic results are shown 
in blue. We see that the agreement is quite good and it is 
better for spheroids with larger elongations. This was expected, 
because the fact that spheroid is strongly elongated was used, 
when deriving our asymptotic formulas, in all the derivations 
and the terms that are small for elongated spheroids were 
neglected. 

For 0 Æ 0, backward wave asymptotics can be found by 
using the same approach, but this has not been done, as yet. 


F. The far field asymptotics 


Since we have been able to obtain a good approximation 
of the induced currents, the use of the Stratton-Chu formula 
in [26] enables us to derive the scattered fields. Since the 
tangential components of the electric vector Æ are zero on the 
surface of perfect conductor, this formula reduces to 


n= -5 ff 3xvaas, 
At 


where J is the total induced current, x denotes the vector 
product and G is the scalar Green’s function 


(56) 


etklr—ro| 


G(r, ro) Pome 

Moving the observation point rp in (56) to infinity along the 
ray defined by the spherical coordinates 79 and Yo we compute 
the limit under the integration sign to obain the formula for 
the far field amplitude of the magnetic field, which reads 


v=- [fax Vyas, 
4r J, 


(57) 


Here 7 is the far field amplitude of G. In view of the 
reciprocity principle, the formula for 7 in the boundary layer 
near the surface coincides with the expression for the field of a 
plane wave incident from the opposite direction. Consequently, 
we need to replace 7 by —7 and y by m — ọ in (34). This 
leads to the asymptotic representations 


+20 
5 Wm cos[m( 
m= 


1 
p= z% + ~ — po), (58) 





gym. e` tkbn+ixn/2 


Viv 


ce ) is 


a 
where Bo = VkbWo is the scaled observation angle, and the 
amplitudes A,, are defined by (36) and (39). 

Our goal is to find the leading order asymptotics of W. 
Towards this end we substitute the asymptotics of the currents 
in (57) from the previous section. 

We use the Cartesian components of the vectors, but perform 
the integration in the coordinates of the boundary layer, where 


dS = b dn ay 1 —- n? dy. 


The expression (57) for the x- and y-componets read 





Vm = 


1 
ttn (59) 
1-7 


Amis, Bo)Mis, oo oe ixT)ds, 


ee Bf (oat) 
v= -4 ff OZ - =F) as. 


The unit vector along the 7 coordinate almost coincides with 
the unit vector along z, that is J, = J, in the leading order, 
but the correction produces the J, component 











Joes eg 
Coe b Tp = n2 n: 
Using the fact that J, = Hy, J, = —H, and the expression 
=i OE, OH. 
H 4+ e) . (60) 
= Sige ie = (27 ? dp 
it is relatively et to get 
a n 
Jz Hy; J, = --—=—H,, 
yp b A — n2 ka 
J = i ( OE, =) 
? kayi- OT Oy 


with the Cartesian components defined by 


Jz = J, cosy — Jy sin p, Jy = J, sin + Jy cos yp. 


Let us now consider the gradient of 7. The differentiation 
by z in the leading order reduces to the multiplication by —ik 
and the other derivatives can be expressed by the formulae 


Ow 1 ((2 Ow 


~ Or 


ðr ayi- 





; ) əy ys 
ixnw } cos p — dp sin Y 
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Fig. 6. Amplitudes of the induced current on the spheroid no. 1 (FEM — red, asymptotics — blue). 


|J] 














0 
-1 


Fig. 7. 


|J] 





-.5 








0 5 


3 























Amplitudes of the induced current on the spheroid no. 2 (FEM — red, asymptotics — blue). 
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Amplitudes of the induced current on the spheroid no. 3 (FEM — red, asymptotics — blue). 
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w z ((2 OV v) si + a cos ) 
= ixn in = ; 
ðy ayl- OT i Op x 
The substitution of the above expressions to the formula for 
W yields 




















1 2T 
v=- fan fap} (252 - OF 2) poos 
Ar Op 
=i, “0 
Ow .. Ow 
Hy (2 siny 4 E cos) b, 
E dE, ôH, 
waz fan fae} (2S De e) ving 
-1 0 





t Hp (23 cos Y a sing) } : 

Next, we compute the integrals in y. The subintegral terms 
contain products of three trigonometric functions. Evidently 
the non-zero contributions are due to terms containing the 
products of three cosines, or two sines and one cosine with 
|m — n| = 1. One can show that in view of the symmetry of 
the problem, the following formulas hold : 











v” = 3 TIF cos(ngo), VFP = X WIE sin(nyo), 
n=1 
wes 3 UTM sin(nyo), vo 2 yM cos(nyo). 


For the harmonics we get 


1 
b O 1 
=g f fyn (2 F 5) (En+1 + En-1) 
=1 


ð 1 
+ (Hanı — Hy-1) (= + 5) vn} dn, 


Vne = 


8 








ð 1 
— (Hn+1 + Hn) (ž + 5) on} dn 


with the plus signs used for the TM incident wave and the 
minus signs to be used for the TE case. 

When substituting the expressions (40) through (43) the 
terms with the indexes n — 2 and n + 2 appear in these 
formulas. With the help of the boundary conditions (27), (28) 
these terms can be excluded. Furthermore, we change the order 
of integration and use the formula 


1 Pew i(t—s) dh 
| (i A 4 1-7? 
=j 
which reduces the integration by s and results in the 
compensation of the terms containing Whittaker functions 





= nd(t — s), 


Myt.n/2(—tx), implying that the incident field gives no conrt- 
ibution. In the other terms the Wronskians of Whittaker func- 
tions (47) are separated, and after substituting the expressions 
given in (36), we get 


+00 
ib(—1)"+1 1 1 
pe et a N 
me m(n!) bBo . 2 2 


x Mit n/2(6B? )M, it np2(i68) (TE 


+00 
i r2 n+l 
2 1 


| + BEE, )dt, (61) 


1 
in)r(* it) x 


x Mitn/2(i8°) Mitn/2 (688) (TEP, — RAP) dt, (62) 
yTM = ib(—1)" + 


+00 
2 n+l Finer n+l , 
TaN i ve ( z 7 u)e( 5 i) x 


x Mit nj2(i8? )M, it wnat iB6) (TM - RIM )dt, (63) 


+00 
ee +1 +1 
on AGRE [eE E-i) 


x Mit nli? )M, it nB (TE 


In the above formulas we have used the fact that TTE = 0 
and T™™ =0. 

We substitute the expressions (44) through (46) and (48), 
for the reflection coefficients in the expressions given in 
(61) throught (63), and rewrite Whittaker functions in terms 
of the Coulomb wave functions. These very cumbersome, 
but straightforward derivations result in the following final 
expressions for the far field amplitudes of the magnetic vector: 


O 
ny = Talbo 

















My RIM )dt. (64) 


+oo 


82b 
op = 3 [fa f do t (bTE H 


+00 


+ X (ae + be + ce + de) cnt) bat, 
{=2 


+00 
8ib E 
v7 BBo f as 





cı + dı) cosy 


(65) 











(66) 


+co 


TM TM 
= —b + cı + d OS Y 
$ TB Bo | f j ns 


+00 


a3 X (a = be + Ce — de) s(t hat 


£=2 


(67) 
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The RCS (in dB) for the spheroid no. | at 1GHz (upper) and 2GHz (lower) for the axial incidence (left) and for the TM wave incident at 3 = 5° 

















Fig. 9. 
(right). 
Fos 1 l \ 1 
Sib de = Fea |t, 58° ) Fess (t, 586) —, 
+00 ee E 1 ; = 1 
- X (a — be — cg + dẹ) cst) bat (68) ye = Hi t 3X Fea | =t, 9X 
l=2 . 1 1 
T —t. —4. 
Here + Hea ( t, 3X Pen ( t DE 
1l 1 1 
ae = Fe-1)/2 (« 5B ) F(e—3) /2 (« 563) = 3 1 , 1 
2 2 2-1 Z= He, (- 5x) Hi (-. 5x) 
1 1 Fo (-t, $x) ; 1 1 
bY = Fo (« ze) Fo (« #8) cae Aid hoe (—t,= 
2 2°°) Hg (-t, $x) ee an ee 
1 1 È E 42 ) We illustrate the approximating properties of the above 
bM = Fy (« 5 g) Fo (« 58) - . ae 7 asymptotics by comparing the RCS 
H =f, aX 
ee RCS = 4n||B||? 
be = Fe- («.56") Poni (« ze) UL. computed by using these formulas and with the help of 
2 2e-1 HFSS of ANSYS!. For this purpose, we choose the same test 
cp = Fe It 1 ii Fea lt 1 B2 Ye Computations were performed by D. Shevnin, CADFEM CIS Branch in 
z z a2 S gO zg North-West Federal Region. 
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Fig. 10. The RCS (in dB) for the spheroid no. 2 at 1GHz (upper) and 2GHz (lower) for the axial incidence (left) and for the TM wave incident at 3 = 5° 


(right). 


problems as specified in Table I. Red and magenta colours 
are used in the Figs. 9 through 11 to identify the results of 
numerical computations and green and blue curves correspond 
to the asymptotic results. 


G. The total scattering cross-section 


The formulas (65) throught (68) provide the asymptotics 
of the far field amplitude of the magnetic field under the 
assumption that the angles V and Wo are small, i.e., the 
directions of the incidence and of the observation are in 
narrow cones near the axis of the body. Nevertheless, these 
formulas enable us to find the total scattering cross-section 
a(o). According to the “optical” theorem [3] we have 


ay = = Im ((p, 90,0), h), (69) 


where h is the magnetic polarization vector of the incident 
wave and the angular brackets denote scalar product. 

The magnetic vector is directed along —e, in the TE 
case and along e, in the TM case. Therefore, we have, 
(PH, htt) =a, Oh) = y7. 
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To exclude the dependence of the effective cross-section on 
the size of the spheroid we normalize it by the visible cross- 


section co = ta? /1 + 82/x and define © = o/oo. Then 


~ nxp? y1 Fax 


+00 +90 +00 
x 5 Im J far + X baa +c + 2 bat (70) 


£=0 





n=0 = 6a 





= x 
Tx B2/1 + B3/x 
+oo +00 +20 
x Xo Im J fima Sin +e zao at (71) 
n=0 6S 


£=0 


For the case of 8 = 0, the representations (70) and (71) 
contain an ambiguity, and by avoiding it we conclude that 
the only nonzero term is co. The formulas then simplify and 
lose the dependence on the polarization of the incident wave; 
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(right). 


hence, we have 


+oo 
1 
fae he w- (12) 
xX zo 1 + e274 


The dependencies of the total scattering cross-section on the 
parameter of elongation x are presented in Fig. 12. We see that 
the rate of elongation may significantly affect the value of the 
scattering cross-section. As y grows, and the body becomes 
less elongated, the effective cross-section tends to its classical 
high-frequency limit of two. 


V. CONCLUSION 


For elongated bodies, the induced currents of the forward 
wave generate the far field amplitudes in the forward cone 
near the axis. For larger observation angles 1), it is necessary 
to consider the currents near the ends of the body. However 
for 3) ~ m the main contribution to Y again comes from 
the currents in the middle part of the spheroids, however the 
backward wave plays a role in this case. The other contribution 
to the back-scattered field is due to the specular reflection. The 
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The RCS (in dB) for the spheroid no. 3 at 1GHz (upper) and 2GHz (lower) for the axial incidence (left) and for the TM wave incident at 3 = 5° 


contribution of the specular point and that of the backward 
wave can both be found in the framework of the approach 
presented in this paper. For the case of acoustic waves, this 
has been done in [27]. 


In addition to spheroids, some other types of elongated 
bodies can also be considered. The parabolic equations enable 
us to use the separation of variables for the case in which 
the quantity p;,(s)p~'/3(s) is constant. This is the case for 
spheroids, hyperboloids, paraboloid and surfaces that have 
no curvature: cones and cylinders that have been analyzed 
in [28]. The case of a paraboloid does not give any new 
results compared to those of Fock [2]. For the one-sheeted 
hyperboloid, some results have been presented in [29]. The 
case of a narrow cone has been studied in [30]. Besides 
the bodies of revolution one can consider infinite cylinders 
with a strongly elongated cross-section. For this case, the 
electromagnetic problems reduce to scalar diffraction and the 
results of [31] are applicable. 

We have considered only the case of plane wave incidence, 
however other sources are possible. For the case of an acoustic 
point source the reader is referred to [32]. The other possibility 
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Fig. 12. The total scattering cross-section as a function of the elongation 
parameter x. Solid black line is for 8 = 0, dash-dotted — for 8 = 0.5, 
dashed — for 8 = 1, dotted — for 8 = 2. Red is for the TE case and blue 
is for the TM case. 





is to consider Gaussian beams incident on an elongated body. 

All the above listed possible directions of further research 
are almost straightforward although they do require derivations 
that are rather cumbersome. 

The other, and more difficult questions pertain to the follow- 
ing two aspects. The first of these is the generalization of this 
technique to the case of the impedance boundary condition, 
making it possible for it to handle dielectric bodies. The 
difficulty here comes from the use the 7 coordinate which 
is related to the normal n by the formula 


ð 1 ð 
OT E 2ay/1 — n? On’ 


which contains 7. The second challenge is to apply the 
approach to bodies whose shapes are not close to one of the 
canonical surfaces, but are composed of parts each of which 
can be approximated by a canonical surface. For example, the 
ogival geometry can be considered as a combination of circular 
cones at the ends and a spheroid in the middle part. We note 
that while asymptotic representations for the diffracted fields 
by cones and spheroids are known, the issue of their matching 
remains an open question. 





REFERENCES 


[1] V. A. Fock “The distribution of currents induced by a pane wave on the 
surface of a conductor”, Journ. of Phys. of the USSR, vol. 10, no. 2, 
pp. 130-136, 1946. 

V. A. Fock “Theory of diffraction from the paraboloid of revolution”, 
in Diffraction of electromagnetic waves by some bodies of revolution, 
Moscow, Sovetskoe Radio, pp. 5—56, in Russian, 1957. 

see also V. A. Fock Electromagnetic Diffraction and Propagation Prob- 
lems, (International Series of Monographs on Electromagnetic Waves), 
Frankfurt, Pergamon Press, Chapter 3, 1965. 

H. Honl, A. W. Maue, K. Westpfahl Theorie der Beugung, Springer- 
Verlag, Berlin—Gottingen—Heidelberg, 1961. 

R. W. P. King, T. T. Wu The scattering and diffraction of waves, Harward 
Univ. Press. Cambridge, Massachusetts, 1959. 


[2] 


[3] 
[4] 


16 


[5] M.G. Belkina, L.A. Wainstein “The characteristics of radiation of 
spherical surface antennas”, in Diffraction of electromagnetic waves by 
some bodies of revolution, Moscow, Sovetskoe Radio, pp. 57-125, in 
Russian, 1957. 

S. Hong “Asymptotic theory of electromagnetic and acoustic diffraction 
by smooth convex surfaces of variable curvature”, Journal of Mathe- 
matical Physics, vol. 8, no. 6, p. 1223-1232, 1967. 

I. V. Andronov, D. Bouche “Computation of the second order term for 
the propagation parameter of creeping waves by boundary layer method” 
Annales des Télécommunications, vol. 49, no. 3—4, pp. 193-210, 1994. 
V. I. Ivanov “Computation of corrections to the Fock asymptotics for 
the wave field near a circular cylinder and a sphere”, J. of Soviet 
Mathematics, vol. 20, no. 1, pp. 1812-1817, 1982. 

I. V. Andronov, D. P. Bouche, M. Duruflé “High-frequency diffraction 
of plane electromagnetic wave by an elongated spheroid”, JEEE Trans- 
actions on Antennas and Propag., vol. 60, no. 11, pp. 5286-5295, 2012. 
T. B. A. Senior “Disk scattering at edge-on incidence”, JEEE Trans 
Antennas Propag, vol. 17, no. 6, pp. 751-756, 1969. 

T. B. A. Senior “Loop excitation of travelling waves”, Can J Phys, 
vol. 40, pp. 1736, 1962. 

T. S. Bird “Comparison of Asymptotic Solutions for the Surface Field 
Excited by a Magnetic Dipole on a Cylinder”, JEEE Transactions on 
Antennas and Propagation, vol. AP-32, no. 11, pp. 1237-1244, 1984. 

I. V. Andronov, D. Bouche “Asymptotic of creeping waves on a strongly 
prolate body, Ann. Télécommun., vol. 49, no. 3-4, pp. 205-210, 1994. 
F. Molinet, I. V. Andronov, D. Bouche Asymptotic and hybrid methods 
in electromagnetics, IEE, London, 2005. 

J. C. Engineer, J. R. King, R. H. Tew “Diffraction by slender bodies”, 
Eur J Appl Math, vol. 9, no. 2, pp. 129-158, 1998. 

S. Yu. Slavyanov, W. Lay, A. Seeger “The classification of Heun’s 
equation and its special nd confluent cases”, in A. Ronveaux (Ed.) Heun 
equation, Oxford University Press, Oxford, 1995. 

F. A. Molinet “Plane wave diffraction by a strongly elongated object 
illuminated in the paraxial direction”, em Progress in Electromagnetics 
Research B, vol. 6, pp. 135-151, 2008. 

N. Ya. Kirpichnikova, M. M. Popov “The Leontovich-Fock Parabolic 
Equation Method in Problems of Short-Wave Diffraction by Prolate 
Bodies”, Journal of Mathematical Sciences, vol. 194, no. 1, pp. 30—43, 
2012. 

I. V. Andronov “High-frequency Asymptotics for Diffraction by a 
Strongly Elongated Body”, Antennas and Wireless Propagation Letters, 
vol. 8, p. 872, 2009. 

I. V. Andronov “High frequency asymptotics of electromagnetic field on 
a strongly elongated spheroid”, PIERS Online, vol. 5, no. 6, pp. 536- 
540, 2009. 

I. V. Komarov, L. I. Ponomarev, S. Yu. Slavyanov Spheroidal and 
Coulomb Spheroidal Functions, Nauka, Moscow, in Russian, 1976. 

M. Abramowitz, I. Stegun Handbook of mathematical functions, Na- 
tional Bureau of Standards, New York, 1964. 

I. V. Andronov “The currents induced by a high-frequency wave incident 
at a small angle to the axis of strongly elongated spheroid”, Progress in 
Electromagnetics Research M, vol. 28, pp. 273-287, 2013. 

I. J. Thompson, A. R. Barnett “COULCC: A continued-fraction al- 
gorithm for Coulomb functions of complex order with complex argu- 
ments”, Computer Physics Communications, vol. 36, pp. 363-372, 1985. 
I. V. Andronov, D. Bouche “Forward and backward waves in high- 
frequency diffraction by an elongated spheroid”, Progress in Electro- 
magnetics Research B, vol. 29, pp. 209-231, 2011. 

J. A. Stratton, L. J. Chu “Diffraction theory of electromagnetic waves”, 
Phys. Review, vol. 56, no. 1, pp. 99-107, 1939. 

I. V. Andronov “High-frequency acoustic scattering from prolate 
spheroids with high aspect ratio”, Journal of the Acoustical Soc. Am., 
vol. 134, no. 6, pp. 4307-4316, 2013. 

I. V. Andronov “Calculation of Diffraction by Strongly Elongated Bodies 
of Revolution”, Acoustical Physics, vol. 58, no. 1, pp. 22-29, 2012. 

I. V. Andronov “The influence of large transverse curvature on the 
induced currents attenuation”, Antennas and Wireless Propag. Letters, 
vol. 11, p. 121, 2012. 

I. V. Andronov, D. Bouche “Diffraction by a narrow circular cone as by 
a strongly elongated body”, Journal of Mathematical Sciences, vol. 185, 
no. 4, pp. 517-522, 2012. 

I. V. Andronov “Diffraction by elliptic cylinder with strongly elongated 
cross-section”, Acoustical Physics, vol. 60, no. 3, to be published, 2014. 
I. V. Andronov “Diffraction of spherical waves on large strongly 
elongated spheroids”, Acta Acoustica united with Acoustica, vol. 99, 
no. 2, pp. 177-182, 2013. 


[6 


oa 


[7] 


[8] 


[9] 


[10] 
[11] 


[12 


= 


[13] 
[14] 
[15] 


[16] 


[17] 


[18] 


[19] 
20] 


[21] 
[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 





[32] 


